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1. Introduction 

In this paper we show that the Bessel equation can be cast, by means of suitable 
transformations, into a system of two parametric oscillator equations, one for a 
damped oscillator, the other one for an amplified one. We show that the group 
contraction mechanism [|l|] is involved in such a relation of the Bessel equation with 
the dissipation/amplification system. 

The interest in such a representation of the Bessel equation is due to the role played 
by the couple of damped/amplified oscillators in several physical problems @. Since 
Bateman [|J proposed to treat the damped system by doubling its degrees of freedom 
and introducing the companion (time-reversed) amplified system, the interest in the 
damped/amplified oscillator system has been growing |4]] and its relevance in treating 
dissipation at classical and quantum level has been stressed f| - @. The couple of 
damped/amplified oscillators has been fruitfully used to describe inflationary models 
of the Universe fLOf , thermal field theories ||, Chern-Simons gauge theory ||11|| , Bloch 



electrons in metals [II], the dissipative quantum model of brain [O, O, Km, etc. It 



also presents features common to the formalism of two-dimension gravity models [15[ 



In particular, the equivalence, under suitable parametrization, between the spherical 
Bessel equation and the damped/amplified oscillators was firstly recognized in the study 
of expanding geometry models Jl0[ and of ordered domain formation in brain models 



On the other hand, the possibility of using and of exploiting the properties of special 
functions in physical problems is in itself of great interest, since, as remarked by Wigner 
16| , "the role which is common to all special functions is to be matrix elements of 



representations of the simplest Lie groups". There is therefore a strong motivation to 
analyze the relation between the Bessel equation and the damped/amplified parametric 
oscillator system also from the perspective of special function theory and group theory. 
We show that the mechanism of group contraction by which one gets the Euclideans 
groups E(2) and E(3), whose representations are given in terms of planar and spherical 
Bessel functions, respectively, introduces the breakdown in the loop-antiloop symmetry 
around a preferred axis. In turn, this can be read off, in a given re-parametrization, as 
the breakdown or violation of the time-reversal symmetry. 

The analysis in this paper brings us also to consider infinite dimensional loop- 
algebras, such as the Virasoro-like algebra, in relation with the Euclidean algebras e(2) 
and e(3). In view of the relevance of infinite dimensional algebraic structures in many 
physical problems, the preliminary results here presented deserve to be further analyzed, 
which is our plan for a future work. 

The plan of the paper is the following. In Section 2 we consider the Bessel equation 
for the spherical Bessel functions and derive from it the set of two damped/amplified 
parametric oscillators. In Section 3 we present similar derivation for the Bessel equation 
for the planar Bessel functions and we discuss the role played by the group contraction 
in our derivation. In Section 4 some preliminary results on the relation between the 
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Virasoro-like algebra and the e(2) and e(3) algebras are presented. Section 5 is devoted 
to conclusions and further remarks. In the Appendix we show that the Bessel-like 
equation describes the damped parametric oscillator under suitable conditions. 



2. The Bessel equation and dissipation 

In this Section we show that the spherical Bessel equation of order n (also called 
the Bessel equation of fractional order |L7]) can be cast, by means of suitable 



transformations, into a set of two equations representing a couple of damped/amplified 
parametric oscillators. In the following Section we analyze the case of the planar Bessel 
equation and the group structure underlying the relation between Bessel equation and 
damped/amplified oscillators. 

The spherical Bessel equation of order n, whose solutions constitute a complete set 
of (parametric) decaying functions fT?], is: 

V 2 Jn;rm + 2 V J n;n + [rf ~ n ( n + 1)] J n = . (1) 

Here n is an integer or zero number, (n = 0, ±1, ±2, ...) and, as customary, the labels 
"; rj" and "; 7777" denote first and second order derivatives, respectively. As well known, 
the solutions of Eq. ([!]), the so called spherical Bessel functions, can be expressed in 
terms of the first and second kind Bessel functions and their linear combinations (the 
Hankel functions). 

Eq. ([!]) is invariant under the transformation n — > — (n + 1) and J n and J_( n+ i) 
are both solutions of the same equation. We can express this by regarding ([!]) 
as an eigenvalue equation and saying that J n and J_( n+ i) are degenerate solutions 
corresponding to the same eigenvalue n(n + 1) of the operator rj 2 + 1r\ + rj 2 . 

We now perform in the Eq.(|I|) the change of variables : 77 — > 77 = ex with x = e _t//a , 
where e and a are arbitrary parameters and the new variable t may be thought to 
denote, e.g., the time variable. By using w n j = J n ■ (x)~ l , Eq.(|l]) then goes into the 
following one: 



21 + 1 . 

U>n,l Wn,l + 

a 



1(1 + 1) - n(n + 1) + e_ 2 e _ 2t/a 
a 2 a 



w n ,i = 0, (2) 



where w denotes derivative of w with respect to time t. We now remark that making 
the choice 1(1 + 1) — n(n + 1) the degeneracy between the solutions J n and J_r n+1 ) is 
removed: in other words, a partition is induced between the two solution sectors {J n } 
and { J_( n+ i)}, as shown by the fact that now two different sets of equations are obtained, 
one for w n j and the other one for iu_/ n+1 w, respectively, each set being composed by 
two different equations, one for I = — (n + 1) and the other one for I = n. The set for 

u> n ,_( n+ i) H w n t n +\) +[{-) e 1 uv-(n+i) = 0, 

a a 

W n ,n ~ Wn,n +[(-)" e~ 2t ^]w n , n = , (3) 

a a 
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for I — — {n + 1) and for I = n, respectively. Similarly, two equations for u>_( n+1 u are 
obtained: 

W-(n+l),n w_ (n+1)in + (-) e 7 ty_ (fH .i) „ = 0, 

a a 

H ,-(n+l) +[(-) 2 e" 2 ' /a ]w_ (rl+ i) _(„+i) = 0, (4) 

for / = n and Z = — (n+ 1), respectively. Inspection of Eqs. (|3|) and Eqs. (f|) shows that 
the symmetry under the transformation n — > — (n + 1) has been broken. 

As a further step, let us now choose the arbitrary parameters a and e to be n- 
dependent: a — > a n and e — > e n (such a choice means that we use r\ — > r\ n = e n x„ with 
x n = e~^ a " ). Also, we perform our choice in such a way that = L and ^- = cuq do 
not depend on n (and on time). By setting u n = w n _( n+ i), and f n = w njH , we then see 
that Eqs. (||) are nothing else than the couple of equations for the damped/amplified 
parametric oscillators: 

U n +Lu n +UJ n 2 {t)u n = 0, 

v n -Lv n +oj n 2 (t)v n = 0, (5) 
with frequency 

u n (t) = uj e"^TT. (6) 



Eqs. (|j) are sometimes called Hill-type equations |l8 |. We see that oo n (t) approaches 
to the time-independent value ojq for n — > oo: the frequency time-dependence is thus 
"graded" by the order n of the original Bessel equation. L and ujq, which may be 
arbitrarily chosen, are characteristic parameters of the oscillator system. 

We note that our choice of keeping L independent of n implies that a_( n+ i) = — a n . 
Then the transformation n — > — (n + 1) leads to solutions (corresponding to J_(„ +1 )) 
which have frequencies exponentially increasing in time (cf. Eq. (|G])). These solutions 
can be respectively obtained from the ones of Eqs. (^) by time-reversal t — > — t and 
exchanging u with v (which we refer to as "charge conjugation"). In the large n limit 
(u n — > ujq) u n and v n are each the time- reversed of the other one and in that limit the 
two sectors { Jj}, i = n, — (n + 1) are mapped one into the other one. 

We thus finally recognize the core of the relation between the spherical Bessel 
equation and the dissipation/amplification phenomenon: the breakdown of the n — > 
—{n + 1) symmetry of the original spherical Bessel equation corresponds in our 
representation to the breakdown of time-reversal symmetry (the emergence of the 
arrow of time) in the manifold of the solutions (the spherical Bessel functions) 
i = n, — (n + 1). 

In the following Section we will see that such a feature has its root in the structure 
of the E(3) group, whose representations, as well known, can be indeed constructed 
by means of the spherical Bessel functions. We will also see that the breakdown of the 
n — > — (n+1) symmetry can be viewed as the breakdown of the loop-antiloop symmetry 
around a preferred direction in the 4D-space. 
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We observe that the functions w n , n and w n r n +i) are "harmonically conjugate" 
functions in the sense that they may be re-conduced to the single parametric oscillator 

r n +n 2 n (t)r n = , (7) 

where r n (t) is defined through w n _ (n+1) (t) = -^r n (t)e^ i and w njTl (t) = -j=r n (t)e%. Vt n 
is the common frequency: 



n n {t) 




(8) 



Remarkably, the first of Eqs. (W), with n — 1 is commonly used in expanding 



geometry (inflationary) models of the Universe JT0[. In that case L denotes the Hubble 
constant. For a discussion of expanding geometry models in terms of Eqs. (^) see |TO . 

In conclusion, the spherical Bessel equation represents, under suitable transfor- 
mations, a two-fold hierarchy of couples of parametric oscillators with constant damp- 
ing/amplification given by L and with time dependent frequency graded by n and by 
—(n + 1). Transition from one tower to the other one (n — > — (n + 1)) is induced by 
time-reversal t — > — t combined with "charge conjugation" u — > v. 

In the Appendix we show that the Bessel-like equation 

j" J= 0, (9) 

where a = 1 or 2, (3 is an arbitrary real constant and ' indicates the derivative with 
respect to the independent variable rj, may represent parametric oscillators with constant 
damping and with different functional choices for the frequency. Vice-versa, from the 
equation of the damped parametric oscillator, by use of convenient transformations, one 
can always obtain the Bessel-like equation (]|). 

In the next Section we consider the group structure underlying the planar and the 
spherical Bessel equation in order to clarify their connection with the damped/amplified 
parametric oscillator system. 

3. Group contraction and Bessel equations 

It is well known that the representations of the Euclidean groups can be constructed in 
terms of the Bessel functions [0, [l?]]- We are therefore interested in these groups in the 
present Section. Features analyzed in the previous Section will be found to be rooted 
in these groups. 

We will focus our attention on E{2) and on E(3), the Euclidean groups in the plane 
and in the space, respectively. Their relation with the Laplace equations in three and 



four dimensions will be summarized (see [16 for details) in view of the connection with 
Bessel equations. 
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3.1. E{2) and planar Bessel equations 

E(2) is the group of the T(v)R(9) transformations, where T(v) is the translation in 
the plane by the vector v (v = (a, b)) and R{9) is the rotation of the plane around the 
origin by the angle 9. The associated Lie algebra is given in terms of the two translation 
generators P a , Pb and of the rotation generator M: 

[P a ,P b } = 0, [P a ,M} = -P b , [P b ,M}=P a . (10) 

The invariant operator of E{2) is P 2 = P 2 + P b 2 = P + P_=P_ P + , with P ± = P a ± iP b , 
which has non-positive eigenvalue, — p 2 . 

The group transformations on square integrable functions / are represented by the 
action of the operator D(v, 9) (see e.g. fl~6H): 

D(v,0)f(<p) = e^/(0), 

D(0,9)f(<j>) = f{(j>-9) . (11) 

In particular, it is possible to assume as basis functions the complete set of 

normalized eigenfunctions of the rotation subgroup {/«(</>)}, f n = (2tc)~^ i~ n e m ^: 

D(0, 9)f n {4>) = e m9 f n ((f)). The complete representation of E(2) can be expressed in the 
form: 

D(v,9)f n =J2 A &8)rnnfm, 
m 

A(V, 9) mn = ^l)m~n e -im P J m _ n{jn .yn{ P -6)^ (12) 

where (r, (3) are the polar coordinates of v and J m is the planar Bessel function of order 
m. Eqs. ( |12D provide the representation of £"(2) in terms of the Bessel functions J m 

0- 

Now we want to study the P 2 eigenvalue equation. To this aim we consider the 
3D-Laplace equation: 

/ d 2 d 2 d 2 \ 

v2 ^(sf + ii + &l)^ ' (13) 

which refers to an isotropic and homogeneous 3D-space. By choosing, instead of the 
spherical or rectangular coordinates, the cylindrical ones, we have: 

(Hit- + A + + JUL) $ = o (i4) 

ydr 2 rdr r 2 d9 2 dx 2 J 
and we search for solutions of the type: 

d 2 



i>(r,0,x 3 ) = (f(r,9) ■ a(x 3 ), —^a(x 3 )=pa(x 3 ). (15) 

ox 3 



Square integrable eigenfunctions (for p positive) are obtained by selecting, for positive 
x 3 , the solution: a = e~ X3P and for negative x 3 the solution: a = e X3P . 

It is clear that the choice of the cylindrical, instead of the spherical coordinates, 
breaks the symmetry of 3D-spatial rotation group SO (3): indeed, when cylindrical 
coordinates are chosen x 3 is differently treated with respect to the two remaining 
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coordinates and this singles out a privileged axis for rotations. The resulting symmetry 
group is E(2), the group contraction of SO (3). This is manifest in the fact that the 
Laplace equation ( |T4"D reduces to the eigenvalue equation for P 2 (realized in polar 
coordinates), i.e. the 2D-Helmholtz equation: 

( d 2 d d 2 \ 2 
+ ^ \dr 2 rdr r 2 d6 2 J ^ ^ ^ ( 

If we solve this equation by assuming <f(r,9) = f(r) ■ e md , we obtain f{r) = J n (pr), 
being J n (pr) the solution of the planar Bessel equation of order n: 

1 n 2 

where rj = pr. Positive/negative values of n correspond to positive/negative rotations 
(loop/antiloop) around the x 3 axis, i.e., they correspond to different orientations of the 



£3 axis. The related solutions are, therefore, different [17]. This occurrence (and the 
similar one in the case of spherical Bessel functions studied below in subsection 3.2) 
bring us in a natural way to consider topological properties of Bessel functions related 
with loop operators and loop-algebras, on which we comment in Section 4. Notice that, 
in spite of the fact that the solutions are not symmetric under the reversal of the x% axis 
(i.e. under time-reversal in our choice 23 = t) Eq. flUTD is invariant under the n — ► — n 
exchange. The n — > —n invariance of Eq. (O) actually reflects the existence of the two 



sets of the £0(3) independent representations: the D n and the D n 1 fH| . 



Now, by using a procedure similar to the one followed in Section 2 for the spherical 
Bessel equation ([I]), by setting w n> i(rj) = (?l)~ l J n (ri), r] = ee~ t//a , and by introducing 
the mirror parameter I = ±n, the planar Bessel equation ( |T7"D gives the couple of 
damped/amplified harmonic oscillators: 



2n . 

I — —n : w n ,_ n H w n n + 

a 

2n . 

I = n : w ntn w n n + 

a 



€ >2 -2t/a 



a 

£\2 e -2t/a 



w n - n = 0, 

Wn.» = 0, (11 



a 

where the value of / is connected with the choice of negative/positive eigenvalues of the 
rotation generator M. 

Summarizing, we have shown that the breakdown of the rotational symmetry of 
50(3) which leads to its group contraction E{2) introduces a crucial difference (loss 
of loop-antiloop symmetry, indeed) in the double choice of the £3 axis orientation. 
This, in turn, results in the difference between the planar Bessel functions or order 
+n and the ones of order —n, in terms of which the E(2) representations can be 
built. Then, we have shown that the planar Bessel equation ( jT7p can be cast, by 



a convenient reparametrization, into the set of eqs. (|18D for the damped/amplified 
harmonic oscillators: the mirror index ±n of the Bessel functions is thus associated to 
the couple of damped/ amplified harmonic oscillators (it is a time-mirror index). 

It is also worth to recall that the contraction of 50(3) to E{2) can be geometrically 
depicted as the projection of the 5 2 sphere (corresponding to 50(3)) on the plane (e.g. 
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the plane tangent to one of the poles of the sphere). It is then clear that the radius, say 
p, of the sphere acts as a "scale" (one may introduce J x = p ■ P a , J y = p ■ P b , J z = M 
with the Js denoting the generators of the algebra so(3)): the E(2) translations in the 
tangent plane are "good" approximations in the limit p — > oo, namely for distances 
much smaller than p. The physical meaning of this is that the 5*0(3) contraction to 
E(2) manifests itself in local observations. However, in the local observation process the 
23 axis orientation is "locked". Which amounts to the loss of symmetry of the solutions 
under the n — > —n exchange (breakdown of the loop-antiloop symmetry). As a matter 
of fact, specifying the direction of the x 3 axis, i.e. choosing one of the two possible 
forms for a (cf. Eq. (|l|)), produces topologically inequivalent configurations ^Oj (see 
also Section 4). 

We also observe that, identifying £3 with the time t coordinate, p has the dimensions 
of an energy over an action (cf., e.g., the eigenfunctions a = e ±X3P ). 

As a final comment, we note that the so called "spinor" representation of S0(3), 
which involves half-integer values of n (corresponding to represent S0(3) in terms of 
the Pauli matrices o~i,i = 1,2,3), maps the planar Bessel equation into the harmonic 
oscillator equation. In fact, for n = ±| and with w = ^/rj J ± i, Eq.(|T7|) reduces to the 
harmonic oscillator equation with frequency p (r plays the role of time): 

^l+p* w (r) = 0, T =X (19) 
ar z p 

In this connection it has been observed |21| that the oscillator so obtained may be 



thought as a possible classical analogue of a Fermi oscillator based on a 'rotation system'. 
3.2. E(3) and the spherical Bessel equation 

The analysis presented in the previous subsection can be extended to E(3), the Euclidean 
group in the space, which is the group contraction of S0(4). The algebra e(3) has six 
generators Pi and Mi, i = 1,2, 3, corresponding to the translation and to the rotation 
generators, respectively. The commutation relations are: 

[^,^=0, [Mi, Mj] = e ijk M k , [P i ,M j \=e ijh P k \ (20) 

We observe that in the contraction process the SO (3) subgroup generated by the M-s 
is left unchanged. The algebra e(3) has two invariants, P 2 = TiPf and EPj • Mj. 

In the 4D-space the Laplace equation for the function ip = ip{x\, X2, £3, £4) may 
be solved by using the position: ip = ip(r,9,(p) ■ (7(24), (r,9,4> spherical coordinates, 
a = e ±XiP ). By a procedure similar to the one of the previous subsection, £4 may be 
considered to play the role of time t. The resulting equation (corresponding to (|TED) is 
solved by the function (p = Y n ^ m (8, <p) ■ J n (pr) where Y n .m is the spherical harmonics and 
J n is the solution of the spherical Bessel equation: 

2 n{n + 1) 

Jnvnri H Jn:ri [1 n ] Jn 0, (21) 
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As in the previous E(2) case for the planar Bessel functions, the spherical Bessel 
functions depend on the continuous eigenvalue p 2 of P 2 and are labelled by n which 
is related with the discrete eigenvalue n(n + 1) of the rotation operator M 2 . 

The order n (integer or zero) classifies the representations D n of the compact 
subgroup SO (3) of E(3). Actually, as already said for the planar case, the existence of 
two sets of 50(3) independent representations (the D n and the D^^ 1 ) is reflected 
in Eq.(f2lf) through its invariance under the transformation n — > — (n+ 1). For n = 0, 
or n = — 1 equation (|2l| ) reduces to the harmonic oscillator equation with frequency p: 

^^+p 2 Mr) = 0, r = 2, (22) 
ar z p 

where w = r\ Jq or r\ J-\. The harmonic oscillator thus appears to be related to the 
"ground state" (in the D n or D^"^ 1 spectra ) of the Bessel system. Here, of course, the 
so-called "true" representation of 5*0(3), i.e., the one with integer values of n has 
been used. 

Again, the breakdown of the symmetry under the transformation n — > — (n + 1) is 
built in the geometrical structure of the E(3) group: the breakdown of the n — ■> — (n + 1) 
symmetry is nothing but the breakdown of the x± axis reversal symmetry (breakdown 
of the loop-antiloop symmetry), i.e. of time-reversal symmetry when x 4 is considered 
to be the time variable, which brings us back to the analysis of subsection 3.1 (see also 
Section 2). Also in the present case, the 50(4) contraction to E(3) manifests itself in 
local observations and the axis orientation then gets "locked". Which amounts to 
the loss of symmetry of the solutions under the n —>■ — (n + 1) exchange. 

4. Loop algebras as extension of Euclidean algebras 

Infinite-dimensional algebras and in particular the so-called "loop-algebras" are 
particularly interesting and useful algebraic structures. It is well known that they can 
be constructed on some finite-dimensional group as SU(1,1) [fl2] and many applications 
exist which exploit such a feature. 

In view of the strict relation between the Euclidean groups, the Bessel functions 
and the dissipation/amplification processes we have discussed in the previous Sections, 
it is also much interesting to consider some topological properties of Bessel functions in 
connection with loop-algebras. 

It is known [|20| that the Bessel functions describe solutions with different 
Pontryagin number in the punctured plane R 2 / (0), where the elements of the homotopy 
group, Il n , may be represented by differential operators acting on analytic functions: 

n n = — , neN, (23) 

with Il„ ■ U m = Il n+m and n is the loop number around the hole. It is possible to 
distinguish two different kinds of behavior, corresponding to two different functions: 

<p m ( Z ) = R m Vm + n{z) , Vn { Z ) = ^ (24) 
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and 

l() m (z) = 4>m-n(z) , 4>m(z) = Z m J m (z) , (25) 

so that on tp, IT n acts in counter-clockwise way while on ip it acts in clockwise way. 
J m (z) is the planar Bessel function (Bessel function of integer order). 

Eqs. (f24|) and ( pop are the well known differential formulae for the planar Bessel 
functions JT7); analogous formulae are true for the functions 

<Pm{z) = Jm(z)z- m , i> m {z) = Z^ j m {z) , (26) 



where the j m are the spherical Bessel functions |20 . 

Notice that while for the planar Bessel functions, the "raising" and "lowering" 
functions coincide for m = (no loops), this is not true with the spherical Bessel 
functions. 

These topological properties, in the planar and in the spherical case, of the ip and 
ip functions readily remind us of the x 3 - and ^-reversal symmetry breakdown discussed 
in Section 3. 

Thus, in view of the topological properties of the Bessel functions and, on the other 
hand, of their relation with the Euclidean groups, it is remarkable that the algebra e(3) 
can be related to a particular structure of loop-algebras. This goes as follows. 

Let us focus our attention in particular on the Virasoro algebra which plays a central 
role in the conformal field theories. We will show that in analogy with the rotation 
algebras, it is possible to follow a contraction procedure which maps the Virasoro algebra 
into a sort of generalization of the Euclidean algebra e(3). This result derives from a 
general construction based on the so called graded contraction method [^. . 

The commutation relations of the Virasoro algebra £ of central charge c (c 
commuting with all the T's) are 

c 

[T n ,T m ] = (n-m)T n+m + — (n 3 - n)5 n+mfl , m,n e Z . (27) 

The Z2-grading of the algebra consists in dividing the set of the T n generators into an 
even set L = {A n , c} and an odd set L\ = {B n }, with 

A n = - ^T 2n + g^n,o^ ) B n = -T 2n+ i , (28) 
so that C = L L\ and 

[La, L ] C L , [Lq, L x ] C L x , [L x , Li] C L . (29) 

The commutation relations of the graded generators are explicitly given by p5| 



2c 

[A n , Am] = (n - m)A n+m + —(n 3 - n)S n+m>0 , (30) 

[B n , B m ] = (n - m)A n+m+ i + ^|(n - -)(n + -){n + ^)5 n+m+lfi , (31) 

[An, B m ] = (n-m - -)B n+m . (32) 
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Eq. (pOD shows that {A n ,c} is again a Virasoro algebra but with central charge 2c. 



We can then consider the Z 2 -graded contraction of the algebra ( pCf ) -fl35|) p5| : 

2c 

[An, A m ) = (n - m)A n+m + —(n 3 - n)5 n+mfi , (33) 

[B n , B m ] =0, (34) 

[A n , B rn ] = (n-m - ^)B n+m . (35) 

Our remark is now that, in the centerless case (c = 0), the Aq and A±\ generators 
close the algebra isomorphic to so(3) ~ su(2) and that the set of these three generators 
and the operators B_i, Bi and B_3 close the e(3) isomorphic algebra. This is shown 
by setting: 

M+= Ai, M - A i , M 3 = iA , 

P + = Bi, P^ = B_i , P 3 = iB_i , (36) 

where the Ms and Ps generators satisfy the commutation relations (p0|) . 

This result has ageneral extension, i.e. the algebra£ n = {A , A± n } {B_i, B ±n _i} 
reproduces the e(3) algebra for each integer value of n, provided the following positions 
are assumed: 

M+ EE -A n , M_ EE -A_ n , M 3 EE -Aq , 

n n n 



P + ee B n _,_ , P_ ee B_ n _ h , P 3 ee iB_ h . (37) 

As final remark we notice that the e(2)-algebra can be obtained as a subalgebra of 
(j3?j) by choosing A± n = 0, for non-zero values of n. 

Our conclusion is that the extension of the Virasoro algebra by means of its Z 2 - 
grading with the subsequent step of the Z 2 -graded contraction appears as a n-graded 
hierarchy of Euclidean algebras. In view of the results obtained in the previous Sections, 
we see that an interesting relation emerges between the couple of damped/amplified 
parametric oscillators graded by n and the loop algebras considered in the present 
Section. 

Exhibiting the link between these mathematical structures appears in itself 
interesting. However, work is still needed in order to fully recognize its physical meaning. 
Further study in such a direction is in our plans. 



5. Final remarks and conclusions 



In this paper we have discussed the relation between the (spherical as well as the 
planar) Bessel equation and the dissipation/amplification processes. Specifically, we 
have shown that, after suitable re-parametrization, the breakdown of the n — > — (n + 1) 
(n — > —n) symmetry of the original spherical (planar) Bessel equation corresponds to 
the breakdown of time-reversal symmetry (the emergence of the arrow of time) in the 
manifold of the solution (the Bessel functions) {Ji}, i = n, — (n+1) {i = n, —n). We have 
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also shown that the mathematical structure through which such a breakdown is realized 
is the mechanism of group contraction leading to the Euclidean groups E(3) and E(2), 
whose representations, as well known, can be constructed by means of the spherical 
and the planar Bessel functions, respectively. For the spherical case, the n — > — (n + 1) 
symmetry of the Bessel equation reflects the degeneracy, for each n, between the two 
independent representations of 50(3), D n and D~ n ~ l . The role of the contraction 
mechanism is the one of removing such a degeneracy. Then, the equivalence between 
the two representations D n and D~ n ~ l is broken and different parametric oscillators 
correspond to each of them. In the sector {i = n} ({z = — (n + 1)},) solutions with 
I = — (n + 1) and with I = n correspond to damped or amplified (amplified or damped) 
solutions, respectively. In the large n limit (u n — > u)q) these solutions are each the 
time-reversed of the other one and in that limit the two sectors { Ji}, i = n, — (n + 1) are 
mapped one into the other. We have seen that, mutata mutandis, similar considerations 
also hold for the planar case. The bottom of the representations (n = 0, or n = — 1 
and in the spinor representation n = ±|) corresponds to the undamped/un-amplified 
harmonic oscillator. 

Summarizing, we can rephrase our finding in the following way: the effect of the 
group contraction (removal of the D n - D^ 11 ^ 1 degeneracy) is to "split" the original 
Bessel equation, conveniently reparametrized, into a couple of damped/amplified 
oscillators for each n. 

We also observe that the fact that the symmetry of the Bessel equation is not 
the symmetry of the solutions reminds us of the phenomenon of the spontaneous 
breakdown of symmetry in Quantum Field Theory (QFT) which occurs when the 
continuous symmetry of the dynamical equations is not the symmetry of the physical 
vacuum. In the Bessel equation case the discrete time-reversal symmetry is broken. 
However, since time-reversal symmetry breakdown manifests itself in the appearing 
of dissipation/amplification phenomena, also continuous time translational symmetry 
is broken (energy non-conservation for the damped (or amplified) system), making the 
resemblance with the QFT case stronger. One further point of contact with the quantum 
case is in the fact that in QFT the effects of spontaneous breakdown of symmetry are 
solely observable in the contraction limit, namely at the observation scale ("locality") 
which is always "small" with respect to the system volume (the infinite volume limit) 
26[ . Also in the present case time-reversal symmetry is broken at local, observational 



level and there dissipation/amplification phenomena become observable phenomena. 

It is interesting to observe that, as well known [Tt\, the Weyl-Heisenberg algebra of 
the quantization procedure is also obtained through the group contraction of the 50(3) 
group (or, in other instances p7[ , of the SU(1, 1) group). 

In this paper our analysis has been limited to the classical sector. The problem 
of quantization of the damped harmonic oscillator, a prototype of dissipative system, 
has been extensively studied in a number of papers B], [[J - |TT[, |0. We recall 
that in order to perform the canonical quantization of the simple damped harmonic 
oscillator, one "doubles" the system by introducing the time-reversed copy of it (the 
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companion amplified oscillator) 0, [|. Canonical quantization is in fact only possible 
for Hamiltonian systems, and the doubling of the dissipative (non-hamiltonian, open) 
system amounts to "closing" it by inclusion of the environment (represented by the 
doubled, amplified oscillator, indeed). At a classical level, however, it is possible to 
solve the damped oscillator equation and, in principle, to treat the dissipative system 
by ignoring the environment (i.e. there is no necessity of introducing the amplified 
oscillator as instead required in the canonical quantization procedure). Remarkably, 
from our discussion in the present paper it emerges that Bessel equation may represent 
both the damped and the amplified oscillator as an inseparable " double-face" unity also 
at the classical level. 

Finally, we have also considered some topological properties of the Bessel functions 
in connection with Virasoro-like loop-algebras and we have presented some preliminary 
results on the relation between such algebras and the algebras e(2) and e(3). More work 
is still necessary in this direction in view of the physical relevance of infinite dimensional 
algebras. 
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Appendix A. Bessel equation vs damped parametric oscillators 



We show in the following that, under several choices of suitable transformations, the 
Bessel-like equation 



t" a t' 

J + - J 

V 



V 2 ) 







(A.l) 



describes damped parametric oscillators with constant damping. In Eq. flA.l ) a = 1 
or 2, (3 is an arbitrary real constant and ' indicates the derivative with respect to the 
independent variable r\. 
By using the position: 

j(n) =x{t) f{ v ), t = t( v ), 

Eq. flA.lp turns into: 

7" «/' n P 2 ] i ' 

x - I - I- 1 — — 

If vf v 2 \(t') 2 , 

Eq.( [A.3| ) describes a damped parametric oscillator, with constant damping 7 and 
frequency u, provided the following two conditions are satisfied: 



fit') 2 {x + x 



f'l t" a 

2- 1 1 

ft>^ {t'Y + T)t> 



0. 



(A.2) 



(A.3) 



2 - — + 



a 



r)t' 1 7 ' 



7 constant, 



(A.4) 
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i.e. 



and 



f" «/' 1 (3 2 ' 



f rjf 
Let us consider conditions 
A) t = 7) , I integer 



rf 



u 2 . 



(|A.6|) for specific choices for t(r)): 



/ oc exp(— ) 77 2 + 2 . 
In the simple case I = 1 (t = rj), u is given by 



1) • 



Note that for ^ — ^- = /3 2 , u 2 gets constant values. 



B) t = exp (qrj), q integer 



cxp 



2 + 2 



with 



7 



o\ 



o + hi- [;(i-o)-^] + 3 



that is, for a = 1 and (3 = \: 



CO 



+ 



gtr? 



1 1 

4 g 2 



1 1 

4 g 2 



1 

f 2 ' 



C) t = In (977), g integer 
/ oc r/^ +1 - Q )/ 2 , 



and 



D) 



lu 2 oc A + Be 2t , ^4, £> constant, 
t = sin (577), g integer 

7 sin gr? _ o/2 1 
/ cx exp — - — 77 7 (cos gr/j 2 ; 



and, for a = 1, (3 = |, 
tu 2 oc .4 



5 



(A.5) 
(A.6) 



(A.7) 



(Ai 



(A.9) 
(A.10) 

(A.ll) 

(A.12) 
(A.13) 

(A.14) 

(A.15) 

ycos qi]) 2 (cosgr?) 4 

It is straightforward that (almost) all the forms for uj are possible. Vice-versa, under 
a suitable transformation, it is always possible to turn the equation for a parametric 
oscillator, with constant damping 7 and frequency u, into a Bessel-like equation of type 
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